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Abstract 

We present a theory for rigorous quantum scattering calculations of probabilities for chemical 
reactions of atoms with diatomic molecules in the presence of an external electric field. The approach 
is based on the fully uncoupled basis set representation of the total wave function in the space- 
fixed coordinate frame, the Fock-Delves hyperspherical coordinates and adiabatic partitioning of 
the total Hamiltonian of the reactive system. The adiabatic channel wave functions are expanded 
in basis sets of hyperangular functions corresponding to different reaction arrangements and the 
interactions with external fields are included in each chemical arrangement separately. We apply 
the theory to examine the effects of electric fields on the chemical reactions of LiF molecules with 
H atoms and HF molecules with Li atoms at low temperatures and show that electric fields may 
enhance the probability of chemical reactions and modify reactive scattering resonances by coupling 
the rotational states of the reactants. Our preliminary results suggest that chemical reactions of 
polar molecules at temperatures below 1 K can be selectively manipulated with dc electric fields 
and microwave laser radiation. 



I. INTRODUCTION 



An important goal of modern chemical physics is to achieve external control over dy- 
namics of elementary chemical processes 1-7 . Manipulating chemical transformations by ex- 
ternal dc fields or laser radiation is at the heart of mode-selective chemistry 1 , chemical 
stereodynamics 3 ' 4 and quantum coherent control of molecular dynamics 6 . External electro- 
magnetic fields can be used to orient and align molecules, which restricts the symmetry of 
the electronic interaction between the reactants in the entrance reaction channel and may 
result in suppression or enhancement of reaction rates, the phenomenon known as the "steric 
effect" 8-10 . Loesch and co-workers 3 ' 11 ' 12 and Friedrich and Herschbach 7 demonstrated that 
rotationally cold polar molecules in the S electronic state can be effectively oriented by dc 
electric fields which was used to study steric effects in molecular spectroscopy 13 , inelastic 
scattering 7 , and chemical reaction dynamics 3 ' 12 . Loesch and Stienkemeier used a combina- 
tion of dc electric fields and infrared radiation pumping to explore the effects of molecular 
alignment in the Li + HF(t> = 1, j = 1) chemical reaction. Their results indicated that side- 
on collisions between HF molecules and Li atoms are more likely to result in the reaction 
than end-on collisions 12 . The steric effects observed in experiments with thermal molecular 
beams are, however, usually weak 3 because the kinetic energy of the reactants greatly ex- 
ceeds the perturbations induced by dc electric fields, even for very polar and heavy molecules 
like IC1 12 . 

Friedrich and Herschbach have shown that molecules can also be aligned by laser 
radiation 14 . The laser alignment method can be applied to both polar and non-polar 
molecules 15 . Stapelfeldt and co-workers 16 demonstrated that significant alignment can 
be achieved with laser fields of 10 12 W/cm 2 . The degree of alignment can be quanti- 
fied by photoionizing the aligned molecules and examining the angular distributions of the 
photofragments 15 ' 16 . Laser-field alignment has been used to produce high-order harmonics 
with specific polarization emitted by N 2 , 2 , and C0 2 molecules 17 . Laser-field alignment 
can also be used to manipulate the rotational motion of molecules 18 or control the branching 
ratios of the photodissociation products 19 . The interaction of molecules with an off-resonant 
laser light is proportional to the square of the electric field strength and substantial alignment 
can be achieved only with very powerful lasers. Because most lasers have short duty cycles, 
laser-aligned molecules are normally produced with low densities insufficient for scattering 
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experiments 4 ' 5 . Other methods, such as collisional alignment in supersonic expansions 4 , pro- 
duce large quantities of aligned molecules, but the degree of alignment in these experiments 20 
is often insignificant and difficult to quantify. 

The effects of external fields on molecular collisions are significantly more pronounced 
at low temperatures. The development of experimental techniques for cooling molecules to 
temperatures near or below 1 K has opened up new possibilities to study controlled chemical 
reactions 21 . Chemical reactions of molecules at cold and ultracold temperatures are accel- 
erated by resonances 22 ' 23 , tunneling 2 ' 22 ' 23 , threshold phenomena 24 , quantum interference 25 
and many-body dynamics 25 ' 26 . With the development of novel experimental methods for 
manipulating molecules with electromagnetic fields such as Stark deceleration 27 ' 28 , magnetic 
or electrostatic guiding 29 and the design of a molecular synchrotron 30 , it has become pos- 
sible to study cold chemical reactions in the presence of external fields experimentally 2 ' 31 . 
Staanum et al. 32 and Zahzam et al. 33 have recently reported measurements of inelastic colli- 
sions and chemical reactions in an optically trapped mixture of Cs 2 molecules and Cs atoms. 
Several research groups are currently developing experiments to study chemical reactions of 
formaldehyde with OH radicals 34 , and Li atoms with HF molecules 35 in slow molecular 
beams. In order to interpret the experimental data and identify new directions for research 
with cold molecules, it is necessary to develop rigorous scattering theory of low-temperature 
chemical reactions in the presence of external fields. 

The first theoretical studies of chemical reactions in the presence of external dc fields 
date back to the work of Karplus and Goldfrey 36 , who used quasiclassical trajectory (QCT) 
calculations based on two different potential energy surfaces to explain the experimental 
observations for the Rb + CH 3 I reaction 8 . More recently, Aoiz and co-workers 37 used QCT 
simulations to explore the effects of electric fields on the reaction DC1 + H — > HC1 + D at 
thermal energies. Aldegunde et al. 38,39 proposed to describe alignment of the reactants in 
terms of polarization moments to examine steric effects in the chemical reactions of H with D 2 
and F with H 2 . Using time-independent reactive scattering calculations, they demonstrated 
that the differential scattering cross sections for the F + H 2 chemical reaction at ultralow 
temperatures can be controlled by changing the polarization of the H 2 molecules 39 . Aquilanti 
and co-workers 40 developed a stereodirected representation for the scattering S'-matrix to 
analyze steric effects in the Li + HF reaction in the absence of external fields. 

Model theoretical studies of chemical reactions in laser fields have been reported by many 



3 



groups 7 . Orel and Miller found that the collinear chemical reactions of H, F, and CI 
atoms with H 2 molecules can be enhanced by intense off-resonant laser fields. Altenberger- 
Siczek and Light 42 used the Floquet formalism to study the collinear reaction LiF + H 
— > LiF + H in an optical field. Orel and Miller 41 reported QCT calculations for the same 
system. Both studies indicated that the reaction probability may be enhanced in the presence 
of laser fields. The enhancement is a result of a laser-induced avoided crossing between the 
ground and the first excited electronic states, reducing the reaction barrier. Seideman and 
Shapiro 44 developed an approximate theory of laser catalysis for chemical reactions in 3D. 
They demonstrated that the chemical reactions H + H 2 44 and D + H 2 45 can be controlled by 
coupling two electronic states with a picosecond laser pulse. Li et a/. 46 have recently applied 
a time-dependent wave packet method to study laser catalysis of the Li 2 + Li exchange 
reaction at low temperatures. Tannor and Rice 47 proposed a mechanism to selectively control 
chemical reactions with an optimized sequence of two femtosecond laser pulses. All of 
the above studies, however, relied on significant approximations to simplify the reactive 
scattering problem in the presence of external fields. For example, the effects of symmetry 
breaking in strong laser fields were not properly considered. These approximations may not 
be adequate for dynamics of chemical reactions at low temperatures. As the collision energy 
of molecules in the subKelvin temperature regime is usually smaller than the perturbations 
due to external fields, the scattering theory of cold chemical reactions must explicitly include 
the interactions with external fields 2 . 

The purpose of the present paper is to develop a rigorous quantum theory of abstraction 
chemical reactions in the presence of an electric field. Our formalism is based on the fully 
uncoupled space-fixed (SF) basis representation of the total wave function expressed in the 
Fock-Delves hyperspherical coordinates to describe the reactive scattering problem. Rigor- 
ous quantum theory of non-reactive molecular collisions in external electric and magnetic 
fields was initially developed by Volpi and Bohn 48 and Krems and Dalgarno 49-51 . These au- 
thors demonstrated that inelastic collisions of molecules at low energies can be manipulated 
with electromagnetic fields. Here, we generalize the work of Krems and Dalgarno to describe 
chemical reactions in external fields. Our calculations indicate that the probability of the 
LiF + H — > HF + Li chemical reaction is sensitive to electric fields of less than 200 kV/cm. 
The remainder of this paper is organized as follows. Section II presents the theory. In 
Sec. Ill, we examine the effects of electric fields on the cross sections and rate constants for 



4 



the Li + EF(v = 0,j = 0) -> LiF + H and LiF(v = l,j = 0) + H -> HF + Li chemical 
reactions at low temperatures. Sec. IV summarizes possible applications of our theory and 
outlines future prospects. 

II. THEORY 

A. Hamiltonian and coordinates 

The Hamiltonian for an atom-diatomic molecule collision system 52 ' 53 in the presence of a 
homogeneous dc electric field can be written in atomic units as 

-[d(R (r a )} • E + H^, (1) 

where R a and r Q are the mass-scaled Jacobi coordinates 53 , fx is the three-body reduced mass, 
7 Q is the angle between the vectors R a and r a , l^(i? a ) r a? 7a) is the interaction potential, 
and £ a is the orbital angular momentum describing the mechanical rotation of the reactive 
complex. In Eq. (1), d(R a , r a , 7 a ) is the dipole moment of the triatomic system, and E is 
the electric field vector which defines the space-fixed quantization axis z. We note that 

lim d(R a , r a , 7«) = d a (r a ), (2) 

where d a (r a ) is the permanent dipole moment of the diatomic molecule. The subscript a in 
Eqs. (1) and (2) refers to different chemical arrangements 52 ' 53 . The asymptotic Hamiltonian 
describes an isolated : E diatomic molecule in the presence of an electric field 50 ' 54 

-^as = H m0 \ + H e f (3) 



where 



and 



1 d 2 i 2 

#moi = — -7r^r a + + V(r a ), (4) 

2fir a dr 2 a 2/ir 2 



H e{ = -d a (r a ) ■ E. (5) 

Here, V(r a ) is the potential energy function, j a is the rotational angular momentum of the 
diatomic molecule, and if e f describes the interaction of the molecule with the electric field. 
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The reactive scattering problem is most conveniently formulated in the Fock-Delves (FD) 
hyperspherical coordinates - the hyperradius p, the hyperangle 9 a , and the Jacobi angle 7 a 
- related to the mass-scaled Jacobi coordinates as follows 



r a = psin6 a , 
R a = pcos9 a . 

The Hamiltonian (1) can be expressed in the FD coordinates 52 ' 53 

£ 2 



(6) 



H = 



1 d 5 d 



+ 



2pp 5 dp dp 2pp 2 cos 2 9 a 



+ V(p, 6 a , la ) - V a (r a ) 

- [d( P , e a , la ) - d a ( P , e a )\ ■ e + h & 



(7) 



The FD asymptotic Hamiltonian can be represented as 



H^b — H mo i + H e {, 



(8) 



where the operators 

-f^mol — 

and 



2/i P 2 







sin 2 26 a d6 a 



sin 2 26. 



d 



+ 



a 89 a sin 2 6 a 



+ V a {r a ) 



(9) 



H e i=-d a (p,9 a )-E (10) 

are obtained by the coordinate transformation from Eqs. (4) and (5). We use the tilde to 
denote the operators and functions expressed in the FD coordinates. 

The total wave function of the reactive complex can be expanded as 52 ' 53 ' 55 



* = p- 5/2 J>(p)<f>^;p), 



(11) 



where ui denotes collectively the hyperangles and p is the hyperradius defined by Eq. (6). 
We emphasize that Eq. (11) is written for a particular component of the wave function 
^i , where i labels the molecular states before the collision. It is therefore important to 
remember that the expansion coefficients in Eq. (11) depend on two indices % and iq, which 
is often written as F iio or F-°. Each column of the matrix of the coefficients F is a linearly 
independent solution of the Schrodinger equation corresponding to the initial state i 52 ' 53 



H<\> = K + \F 

11 ^ IQ J -'t ^IQI 



(12) 
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where E t is the total energy. We will use the index i Q only when necessary. 

The basis functions $j(cj;p) can be chosen as the eigenfunctions of the adiabatic 
Hamiltonian 53 ' 55 

H &A (p)$i(u; p) = ei(p)<S>i(u; p), (13) 

where €i(p) are the adiabatic eigenvalues, which depend parametrically on p. In this work, 
we define the adiabatic Hamiltonian as the total Hamiltonian (7) without the hyperradial 
kinetic energy 

^ad(p) = HM + H e{ , (14) 

where H e { is the Hamiltonian for the molecule-field interaction given by Eq. (10) and the 
field-independent term can be written in the form 

= o 2n + V(p, 9 a , 7a ) - V a (r a ) + if mol . (15) 

COS U 

In deriving Eqs. (14) and (15), we have neglected the interaction of the three-body compo- 
nent of the dipole moment with electric fields 

#ef,3B = ~[d(p, 9 a , 7a) " d a (p, 9 a )\ ■ E. (16) 

Ab initio calculations show that the expression in the square brackets decreases rapidly with 
increasing p 56 . At small p, the interaction of a transient reaction complex with electric 
fields is negligible in comparison with the electronic interaction potential of the complex 57 . 
Since the interaction due to the three-body term in Eq. (16) couples the same states as the 
interaction potential, it can be safely neglected 57 . It was demonstrated in our previous work 58 
that including the interaction with the transient dipole moment in scattering calculations 
does not modify the collision dynamics except near certain scattering resonances at extremely 
low temperatures. 

Substituting the expansion (11) into the Schrodinger equation with the Hamiltonian (7) 
and using Eq. (13), we obtain a system of close-coupled (CC) equations 52 ' 53 ' 59-61 

d 2 15 



+ 2p[E-e i (p)} 



Hp) = o. (17) 



dp 2 8pp 2 

Eqs. (11), (13) and (17) can be solved on a grid of radial sectors p k extending from small 
values of p to p = oo to determine the total wave function of the reactive complex subject to 
the scattering boundary conditions. The reactive scattering problem can thus be separated 
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into three steps: (i) the diagonalization of the adiabatic Schrodinger equation (13); (ii) 
the integration of the CC equations (17), and (iii) matching the results of the numerical 
integration to the asymptotic form of the wave function as determined by reactive scattering 
boundary conditions. 

B. The adiabatic eigenvalue problem 

The adiabatic eigenvalue problem (13) can be solved using an expansion 

$^;p) = ^T m |C„>, (18) 

n 

where \( n ) are some orthonormal basis functions. Substituting this expansion into Eq. (13) 
leads to a matrix eigenvalue problem 

^[<Cn|#ad(p)|Cn'> ~ Q(p)<W] T n'i = 0. (19) 
n' 

The eigenvalues e^(p) and eigenvectors T ni can be found once the adiabatic Hamiltonian 
matrix is specified in the basis In order to construct this matrix, we introduce a basis 
set that simultaneously includes primitive functions defined in all chemical arrangements 62 ' 63 

Xavj (0 a ; p) \jMj) \£M e ) = Xavj {0 a ; p)Y ]Mj {r a )Y eMe (R a ) (20) 

where the uncoupled space-fixed (SF) basis functions are direct products of the spherical 
harmonics Yj M .(f a ) and Yi Me (R a ), and the FD rovibrational eigenfunctions and eigenenergies 
are defined as 

H m0 \Xavj(6a', P) = £avj(p)Xavj(9a', P) ■ (21) 

The projection of the total angular momentum onto the electric field axis M = Mj + Mi is 
rigorously conserved for collisions in parallel fields 48 ' 49 . Therefore, the basis functions (20) 
corresponding to different values of M are not coupled, and the Schrodinger equation (14) 
can be solved independently for each M. 

We emphasize that the basis functions (20) are not the eigenfunctions of the total angular 
momentum operator and that the angular momentum projections Mj and Mi are defined 
with respect to the SF quantization axis determined by the direction of the external field. 
All previous methods for solving the reactive scattering problem used the total angular 
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momentum representation of Arthurs and Dalgarno 64 in the body-fixed (BF) coordinate 
system. The quantization axis in the BF approach is directed along the Jacobi vector R Q . 
This choice simplifies the evaluation of the matrix elements of the interaction potential. In 
addition, chemical reactions with near-collinear transition states are determined by a limited 
number of the BF projections of the total angular momentum 65 , which leads to a substantial 
reduction of the number of scattering channels. External electric fields break the isotropy 
of space and couple states corresponding to different total angular momenta and inversion 
parities. Therefore, the total angular momentum representation 64 ' 66 offers no advantage. 
Krems and Dalgarno showed 49 that in the presence of an external field, it is more convenient 
to work directly in the SF frame as this leads to simpler expressions for the matrix elements 
of the molecule-field interaction. Although the matrix elements of the interaction potential 
in the SF basis are more complicated 49 , they can be obtained from the BF matrix elements 
using a simple transformation as shown in Sec. II C below. 

We note that the basis (20) is not orthogonal because the functions Xavj(9 a ', p) of different 
a overlap at small p 53 . An appropriate orthogonal basis set can be defined in terms of the 
functions (20) following the symmetric orthogonalization procedure 65 ' 67 ' 68 

|Cn) = "i= X c m jM j lM t , n Xauj{Qa,p)\jMj)\tM t ), (22) 

V n 3Mj,tMl 

where \ and X^ v - M , m ^ n are, respectively, the real eigenvalues and eigenvectors of the 
overlap matrix 

X T OX = A, (23) 

where A = diag(Ai, . . . , A n ) and the symmetric and orthogonal overlap matrix of the primi- 
tive basis functions (20) is given by 

OavjM j eMi,a'v'j'M' j e'M' e = 

(x^<^;p)IO'M,|(£M,| XaV/ (^;p))|/Mj)|fM;). (24) 

It is easy to verify that the functions defined by Eq. (22) are orthogonal, and therefore 
can be used to expand the adiabatic surface functions in Eq. (13). The matrix elements of 
the adiabatic Hamiltonian (19) can be obtained from Eq. (24) 

(Cnl-^ad(p)ICn') = 7= = ^ ^ ^avjMjiM^n^a'v'j'M'/'M^n' 

V n n ' a,v,j,Mj,t,Mt a 1 ,v' ,j> ,M'j,Z' 

x (x^(^;p)|(jM J -|(£M,|// ad (p)| XaVj v(^,;p))|/M')|/M;). (25) 
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Because the adiabatic functions calculated at different p are not orthogonal, we need 
an additional transformation between the adiabatic functions of the adjacent propagation 
sectors (say, pk-i and pk). The transformation is derived in Sec. HE and has the form 

[$>(Pk-i, Pk)]u> = ($i(w;pfc_i)|$i'(w;p fc )), (26) 

where the integration is carried out over all variables except p. Expanding the adiabatic 
surface functions as in Eqs. (18) and (22), we find 

[S(pfc-1, Pk)]ii' — T ni T n iji - XavjM i lM l ,n{,Pk-^X-cx l v'j l M' i t'M[,n l {,Pk) 

n,n' * n a,v,j,M jt e,M e a'v'j'M'/'M^ 

x (Xa^(^;p fc -i)|(iM J |(£M,| XaV/ (^,; Pfe ))|/A^.)|fM;) (27) 

Unlike the overlap matrix calculated at fixed p (24), the overlap matrix of the primitive 
functions on the right-hand side of Eq. (27) is not symmetric. We denote this matrix as 
Osf(pa;-i, Pk)- It is discussed in more detail in the next section. 

In summary, the adiabatic eigenvalue problem can be solved in three steps. First, the 
overlap matrix of the primitive basis functions (24) is constructed and diagonalized to yield 
the matrices A and X. Second, the matrix of the adiabatic Hamiltonian in the primitive basis 
is evaluated and transformed to the orthogonalized basis as in Eq. (25). Third, the matrix 
(25) is diagonalized to give the eigenvectors T ni (p), which are convoluted with the sector- 
to-sector overlap matrix (27) to yield the sector-to-sector transformation matrix S(pk-i, Pk)- 
In order to implement this strategy, we need to evaluate: (i) the overlap matrix between the 
primitive functions (24); (ii) the matrix of the adiabatic Hamiltonian (25) in the primitive 
basis, and (iii) the sector-to-sector overlap matrix Osp(Pk-i, Pk) in the primitive basis. 

C. Matrix elements in the primitive basis 

1. Overlap matrix 

Because the primitive functions of the same arrangement are orthonormal, we have 

QavjMjtMt&v'j'M'jl'M't = ^vv'^jj'^MjM'^W^MtM^ (28) 
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The only nonzero elements of are those between different arrangements. They have the 
form 

tt/2 



/r r 7T / z 
dR a J dr a J d9 a YjMj{f a )Yji M i.{f a i) 



x Y mt (R a )Y VM >(R a ,)[ . a )<j> av j(O a ; p)<f)a'v'j'(6a'', p), (29) 

y Sin Zu a i J 

where the renormalized FD rovibrational basis functions (f> av j(9 a ; p) are defined by 

x a ^p) = 2KvA l; p) - (so) 

sin zo a 

The functions in the integrand (29) depend on different coordinates. A coordinate trans- 
formation is therefore required to express the functions of one arrangement in terms of the 
gridpoints of the other. The FD coordinates of different arrangements are related by kine- 
matic rotations 52 ' 53 ' 69 . These rotations lead to a complicated inseparable dependence of the 
vectors R a /, f a i on the integration variables R a , f a . The integral in Eq. (29) is thus a 
five-dimensional integral, and is computationally intractable. 

The problem can be simplified by transforming Eq. (29) to the BF frame, which allows us 
to integrate over the three Euler angles analytically. The transformation of the SF functions 
has the following form 52 ' 53 ' 70 

Xavj(O a ;p)Y jM .(f a )Y tMt (R a ) = Xavj (6 a ,p)(-) M (2£+1) 1 / 2 

j+i ( ■ p j \ min(j,J) (ad J \ 

x £ (2J+1)W 31 £ [ 3 )\JMK)\jK), (31) 

where the BF basis functions are given by 

(1 J + 1\ 1/2 

Xavj (0a;p)\JMK)\jK) = Xa Vj (0a;p) ( -gj" ) D J MK(VE)V2^Y jK ( la ,0). (32) 

Here, K is the projection of the total angular momentum J along the BF z-axis defined 
by the vector R a , the symmetric top functions \ JMK) depend on the Euler angles fi E 71 ' 72 , 
and the renormalized spherical harmonics \jK) = \p2^YjK{la, 0) describe the rotation of 
the diatomic molecule in the BF frame 63 . 

Using Eq. (31) and the orthogonality properties of the BF functions, we can rewrite the 
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SF overlap matrix (24) in the form 

O av jM j Me,a'v'j'M' j e'M' t = [{21 + l)(2£' + 1)] 1/2 ^(2J + 1) 



X 




3 t J 

Mj M e -M ) \ M\ M[ —M' J 

x (xavAp^a)\(JMK\(jK\ Xa 'MP^a')\JMK')\j'K'), (33) 

The matrix elements on the right-hand side are much easier to evaluate because the overlap 
between the symmetric top functions is (JMK'\JMK) = dx K ,(uj aa >), where the angle u aa > 
between the vectors R a and R a > is a function of two angles 9 a and 7 Q 53 ' 62 . The result is well 
known 53 ' 62 ' 63 and we present it here for completeness 

/»7r/2 rir 

(Xav j (P^a)\(JMK\(jK\ X a'v' j ip;0 a/ )\JMK')\j'K')^27r / d6 a / sin la d la 

Jo Jo 

x ( a I Y jK(la, O)(f> avj (0 a ; p)d J KK ,(uj aa/ )Y j < KI (>y al , 0)(f> a > v >j>(9 a r, p) (34) 

ysm Zu a i J 

2. Adiabatic Hamiltonian and sector-to-sector overlap matrices 

The matrix elements of the adiabatic Hamiltonian (14) can be derived following the same 
procedure. The field-independent part of the adiabatic Hamiltonian (14) can be written in 
the BF frame as 70 

& M = 9 (J 2 ~ Jt £ + V (Pi 7a) - V a (r a ) + H moh (35) 

ZilXQ COS (/q/ 

This Hamiltonian is formally equivalent to that given by (14). Therefore, we can evaluate 
the matrix elements of H&{p) in the SF frame using the transformation (31) as follows 

(Xavj(O a ; p) I (jM J I (£M e \H 6 (p) \ Xa , v/j , {9 a , ; p)) |/Mj> = [(2£ + 1)(2£' + 1)] 1 / 2 ^(2J+ 1) 

3 t J \ I f V J 



x 



y^{-i) K+K ' 

Mj M e -M J \ M' Mi —M' ' 




x (Xa V M°a)\(JMK\(jK\H 6 (p)\x a 'MP^^\JMK')\j'K'), (36) 

Since the interaction potential does not depend on the Euler angles, the matrix elements of 
the field-independent part of the adiabatic Hamiltonian in the BF basis can be reduced to 
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2D integrals over 9 a and 7 Q . They have the form 53 ' 63 ' 65 

(Xa V M^)\(JMK\(jK\H 6 (p)\ Xa ' v 'f(p;e a/ )\JMK%'K'} = 2tt I ' ' d6 a [ sin 7a d 7 , 

Jo Jo 

3 ~ Jm) " +V(p,9 a , la )-V a (r a ) + H mol 



x Y jK (>y a ,0)(f) avj (6 a ; p)d KK ,( 



2pp 2 cos 2 9 a 



sin 29 a 

X ( • of) a ) Y j'K>(la>,0)<f>a>v>j>(O a >; P)- (37) 
Sill 2aj Q/f 



The matrix elements of the first term in square brackets can be obtained by expressing the 
angular momentum operator in terms of the raising and lowering operators, and taking into 
account the anomalous commutation relations as described in Refs. 52 ' 53 ' 62 ' 71 . The remaining 
terms in Eq. (37) are very similar to the overlap matrix elements discussed in the previous 
section. 

Finally, the sector-to-sector overlap matrix in the primitive basis has the form 

(x^<^;p fc -i)l(jM,|(£M £ | XaW (^,;p,0)|/M;)|£ / M;) = [(2£+ 1)(2£' + 1)] 1 / 2 ^(2J+ 1) 

j 



x 



Mj M e -M J \ Mj M; -M' J \K0-kJ\K'0 -K' 



x (xavApk~i,0«)\(JMK\(jK\xa'v>f(pk;O a ,)\JMK')\j'K'}, (38) 
where the BF sector-to-sector overlap matrix is given by 

/ i 7r/2 /»7T 

(Xavj(Pk-i;O a )\(JMK\(jK\x a ' v 'f (p k ;0 a ,)\ J MK')\fK') = 2tv / dO a / sin 7a rf 7a 

Jo Jo 

X ( -T-^r WjKilaity^ocvjiO^ Pk-l)dK K ,{u aa ')Y? (39) 

D. Interaction with electric fields 

If the interaction of the three-body term with electric fields is neglected, Eq. (10) can be 
rewritten as 

H ei = -d a (p, 6 a )E cos Xa, (40) 

where Xa is the polar SF angle of the diatomic molecule in arrangement a. In the limit of 
large p, there is no overlap between the vibrational basis functions of different arrangements, 
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and the matrix elements of the interaction with electric fields (40) can be written as 



{Xavj(0 a ; p)\(jM j \(eM e \H e t\x a >v>A d °'; P))\f M j)\^ M e) = -E5 aa/ 5 ee 5 MiM ' e 

X (XavjiOa, P)\d a (p, a )\Xav'j'(0 a ; P))(J M j\ COSXalfM'j). (41) 

Furthermore, at large p the dipole moment function d a (9 a ; p) is localized near the equilibrium 
distance of the diatomic molecule and it is a slowly varying function of r a . Therefore, it is 
a good approximation to neglect the matrix elements off diagonal in v: 

(Xavj(0 a ; p)\d a (p,9 a )\x av 'j'(0 a ; p)) = 8 vv >d av . (42) 

We note that this approximation is consistent with neglecting the coupling between different 
arrangements due to the electric field (41). For low v considered in this work, the vibrational 
dependence of the matrix element in Eq. (42) is very weak. For example, the dipole moment 
of LiF in the v = 8 vibrational state is only 11% larger than for the ground vibrational 
state 82 . We therefore assume that the matrix element in Eq. (42) is equal to the dipole 
moment of the molecule in the ground vibrational state and is independent of v. Writing 
cosxa = (47r/3) 1 / 2 Fio(r a ) and evaluating the integral over the product of three spherical 
harmonics 71 ' 72 , we find 48 ' 50 

(Xavj(0 a ; p)\{jM j \{eM e \H e f\xa>v'A d «'l p))\j'M^)\e'M' e ) = -Ed a 5 aa ,5 vv ,5 u ,5 Me M' e 

x(-)«[(2j + l)(2j' + !)?'*( j 1 f ) ( ilf ), (43) 

y —Mj M'j J y J 

where the symbols in parentheses are the 3—j symbols. In order to complete the definition 
of the matrix elements of the adiabatic Hamiltonian (14), we need to evaluate the matrix 
elements of the asymptotic Hamiltonian (8). Because the overlap matrix (33) becomes 
diagonal in the asymptotic region, the matrix elements of the asymptotic Hamiltonian can 
be evaluated using Eq. (43), 

(Xavj(0 a ; P) I {jMj I (£M e \H as | XaV/ (0 a ,; p)) l/Mj) \t'M[) 

J 1 f \ ( j 1 f 



X 



S,A M]M/j e avj (p) - Ed a (-) M *[(2j + l)(2j" + I)} 1 ' 2 



-Mj M'j j \ 



(44) 



This expression shows that the asymptotic Hamiltonian is not diagonal in the uncoupled SF 
basis (20), except for E = 0. The scattering boundary conditions must be applied in the 
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basis that diagonalizes the asymptotic Hamiltonian (44) 49 . The eigenfunctions of H as 

H AS ifjavrM T (r a , a ;p) = e aVT M T {p)ipavrM T (r a , 6 a ; p) (45) 
can be written as linear combinations of the basis functions given by Eq. (20) 

ipavrM T (r a , 9 a ; P) = ^2 C jM J ,rM T (E)xavj(O a ; p)Y jMj (r a ) , (46) 

where the field-dependent mixing coefficients CjM 3 ,TM T (E) can be obtained by numerical 
diagonalization of the asymptotic Hamiltonian matrix (44). Since the electric field only 
couples the levels with Aj = ±1, the dependence on j of the overlap of rovibrational FD 
basis functions in Eq. (44) can be neglected, and Eq. (46) rewritten as 

XavAOa, P) ^ C jM J ,rM r (E)Y jMj (f a ), (47) 

where r, M T labels the field-dressed states. The index r of the function Xavri^a', p) denotes 
the dominant rotational state in the expansion (47). 

The matrix of the asymptotic Hamiltonian (3) can also be written in the Jacobi basis 
defined in the same way as the FD basis (20) 

^M\ jMj )\£M e ) : (48) 

where the functions \jMj) and \£Me) are the same as in Eq. (20) and the renormalized 
Jacobi rovibrational basis functions are the radial eigenfunctions of the diatomic molecule 
(4) in the absence of an electric field 



+ o 2 + V ( r «) 



2p drl 2prl 
The eigenfunctions of H m0 \ in Eq. (3) are given by £, av j(r a )/r a . 

Although the basis (48) is slowly converging at small p (where different reaction arrange- 
ments are strongly coupled 52 ' 53 ' 62 ), it is convenient in the asymptotic region, where the Jacobi 
basis functions of different arrangements are orthonormal. The matrix representation of the 
asymptotic Hamiltonian (10) in this basis is given by 



(Za»j ij a ) | (jMj | {£M e \H as \^ vlj/ (rv)> I J >M;) \£'M' e ) = 5 aal 5 vv i5 U '5 Ml M' t 

j 1 f \ I j 1 f 



x 



S jf d MjM ^ avj - Ed a (-) M i[(2j + l)(2j" + l)] 1 / 2 

J J 



-Mj M'- I \ 



, (50) 
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where e av j is the rovibrational energy of the diatomic molecule defined by Eq. (49). 

Since e av j(p) — > e av j in the limit of large p, Eqs. (50) and (44) define the same matrix, 
and the eigenvectors CjM J>T M T m Eq. (47) are identical in the Jacobi and FD coordinates. 
This is important for reactive scattering boundary conditions (see Sec. II F). 

E. Propagation 

In order to determine the expansion coefficients Fj(p), it is necessary to integrate the CC 
equations (17) from small values of p to the asymptotic region. Since the adiabatic basis 
(13) changes with p, the solutions Fi{p) should be transformed to a new basis as p increases. 
The integration interval of p is usually divided into small sectors such that the adiabatic 
basis (13) does not change within the sector. At the boundary between the (k — l)-th and 
/c-th sectors, the solutions of Eq. (17) must be transformed to the basis of the /c-th sector. 
The transformation is determined by the requirement that the total wave function (11) be 
continuous at the boundary. The coefficients -Fi(Pfc) can be written as 

Fi(pk) = y^Fj/(p fc _i)[S T (p fc -i,Pfc)]n', (51) 

v 

where pk-i and pk denote the centers of the respective sectors, and S(p k -i, p k ) is the sector- 
to-sector overlap matrix given by Eq. (26). 

To preserve the numerical stability of the solutions of the CC equations, we propagate 
the log-derivative matrix 52 ' 53 ' 73 ' 74 defined as 

Y(p k ) = F'(p k )[F(p k )}-\ (52) 

where F(p) is the matrix with columns represented by the solutions to Eq. (17), and the 
prime indicates differentiation with respect to p. We propagate the matrix Y across the 
sector [pk-2, Pk-i] (the first sector of the grid is [po,pi]) using the diagonal reference potential 
propagator of Manolopoulos 74 . After reaching the right end at p k -i, we transform the log- 
derivative matrix to the adiabatic basis of the next sector [p k -i,p k ]. Eq. (51) can be written 
in matrix form as 

F(p k ) = S T (p fc _i,p fe )F(p fe _ 1 ). (53) 
Using this relation and the definition (52), we obtain 

Y(p fc ) = S T (p fe _i,p fc )Y(p fc _ 1 )S(p fe _i,p fc ). (54) 
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The CC equations can thus be integrated by repeated application of the log-derivative 
propagation 74 and the transformation (54). 



F. Asymptotic Boundary Conditions 

In the asymptotic limit p — > oo, different reaction arrangements become uncoupled, and 
the reactive scattering wave function can be conveniently re-expressed in the SF Jacobi basis 

*= E E 

(R a )il>avTM T (r a )YeM e (R a ), (55) 

a,v,r,M T i,Me Ta a 

where the asymptotic field-dressed Jacobi functions of the diatomic molecule are given by 

HaslpavrMri^a) — ^avTM T ^J J avTM T ( r a)- (56) 

The energies e avT M T are the Stark levels of the diatomic molecule in arrangement a. The 
field-dressed eigenfunctions in the Jacobi coordinates defined by Eq. (56) are distinct from 
the functions in the FD coodinates (45), and we will specify the arguments of both functions 
to avoid confusion. Similarly to Eq. (46), the Jacobi asymptotic functions can be expanded 
in spherical harmonics 

i>avTM T (r a ) = (, aVT (r a ) ^ Cj Mj ,TM T {E)Yj Mj (r a ) (57) 

with the expansion coefficients CjM J>T M T {.E) given by Eq. (47). In Eq. (57), £ a vj( r a) is 
the Jacobi ro-vibrational eigenfunction in arrangement a given by Eq. (49). The radial 
functions in Eq. (55) have the following asymptotic behavior 52 ' 66 

F aWM™' e M' e ( R <x' -» oo) ~ S aa i5 vvl 5 TTl 5 MT M^wSM e M' t exp[-i(k avr M T R a - £n/2)] 

( KvrM \ l/2 

— I -j — I S a 'v'T'Mie'MkavTM T eM t e^[i{ka'v'T'MlRa' — ^ 7r /2)] 

(58) 

where k 2 avrMT = 2fx(E t — e aVT M T ) is the asymptotic wave vector given in terms of the total 
energy E t and the asymptotic Stark energy (56), and Sij are the ^-matrix elements. They 
can be evaluated by matching the asymptotic form of the FD wave function (11) to the 
Jacobi wave function (55). The details of the asymptotic matching procedure are described 
in the Appendix. 
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The scattered part of the wave function (55) corresponding to the initial flux of molecules 
in the state (a, v, r, I) propagating along the direction R a . has the form 52 ' 53 ' 66 

^avTM^i^-a' — » OO) ~ ^ i(k a vTM T k a 'v' V M;)" 1 ^ QavrMr^a'v'r' M' T {Ra v Ra') 

a! ,v' ,t' ,M' T 

exp[ik a > v ' T /M'Ra>] , , s ,_„>. 

x Wa'v'r'M' T {r a ), (59) 

r a >K a > 

where q aVT M T ^a'v'T'M^(R ai , Ra') is the the scattering amplitude. An expression for the scat- 
tering amplitude in terms of the ^-matrix elements can be obtained by substituting Eq. 
(58) into Eq. (55). After separating the scattered component of the wave function and 
comparing the result with Eq. (59), we find 49 

QavrMr^a'v'T' M^{Ra v Ra') = 27T ^ ^ ^ i Y e * Me (R a .)YiM e (Ra') 

e,Me I'M' 



X 



\^aa'^vv'^TT'^M T M' T ^lt'^M t M' l ~ S avr M T £M e ;a'v'r' M' T i' M'^\ ■ (60) 



The differential cross section for reactive transitions between the field-dressed states is given 

by 

(R ai ,R a ,)\ 2 . (61) 



dR a . dR a i kavrM T 

Integrating this expression over R a / and averaging over R a . gives the integral reaction cross 
section 

VavTM T ^a'v'T'M' T = p ^aa'^vv'^TT'^MrM^U'^M^ ~ S avr M T £M e ;a'v'T' M' T t' M'^ ■ 

(62) 

Taking into account the conservation of the total angular momentum projection M, we 
can rewrite this expression as a sum of partial cross sections calculated from the S'-matrix 
elements at fixed M 

OavrM T ^a'v'T'M' T = p ^2 \^aa'^vv'^TT'^M T M^U'^M t .M' i ~ ^avrM-rlM^a' v' t' M' T i' M'^ ■ 

avrM T M e,Met',M^ 

(63) 

where the summation over Me and M[ is restricted so that Mi + M T = M[ + M' T = M. Eqs. 
(60)-(63) generalize the expressions derived by Krems and Dalgarno 49 to reactive scattering 
in external fields. In the absence of an external magnetic field, the Stark levels with energies 
tavrMr and e aVT -M T are degenerate. Because the sum j + f + 1 in Eq. (44) is always even, 

18 



the matrix elements of the asymptotic Hamiltonian do not depend on the sign of M and so 
do the M-resolved cross sections given by Eq. (63). With this in mind, Eq. (63) can be 
rewritten as 

CavrMr^a'v'r'M^ = p ^ ^ ^ (2 — S M fl) 

avrM T M >0e,M e £',M' e 

X \S aa 'S vv >S TT iSM T M^UiS Me M' e - Sa VT M T eM e ;a'v'r'M^'M' e \ 2 (64) 

III. RESULTS AND DISCUSSION 

In this section, we present the results of preliminary calculations based on the theory 
described in Sec. II. We study the variation of the cross sections for the chemical reactions 
LiF(w = 1, j = 0) + H -> LiF + H and Li + RF(v = 0, j = 0) -> LiF + H with the collision 
energy and the strength of an applied electric field. We also analyze the competing process 
of vibrational relaxation in LiF (v — — 0)+ H collisions and suggest a mechanism for 
electric field control of chemical reactions at low temperatures. 

A. Computational details 

The hyperspherical basis functions Xavj(0 a 'i p) are constructed by solving Eq. (21) using 
the Fourier grid Hamiltonian method 75 . The functions become more localized with increasing 
p. The localization region of the functions (21) can be defined as follows 65 ' 76 

C in = sin" 1 [C n /p] ; C ax = sin" 1 I P\ ■ (65) 

We chose the endpoints r™ m and r™ ax so that the rovibrational function of the diatomic 
molecule (49) is nonzero in the interval [r™ m ,r™ ax ]. This guarantees that the functions 
Xavj(9 a ',p) vanish outside the interval [#™ m ,0™ ax ], and the integration range in Eqs. (34) 
and (37) reduces from [0, 7r/2] to [8™ m ,6™ ax ]. In addition, using the p-dependent endpoints 
(65) ensures the accuracy at large p which is necessary to calculate reactive scattering cross 
sections at low temperatures. 

The BF matrix elements given by Eqs. (34), (37) and (39) are evaluated at fixed J 
as described by Miller 62 , and Alexander, Manolopoulos and Werner 63 . The integration is 
performed with 45 Gauss-Legendre quadrature points in 9 a e [#™ m (p), #™ ax (p)] and 7 a G 
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[0, 7r]. The basis functions are transformed between different reaction arrangements using 
the expressions derived in Refs. 53 ' 62 . The upper limit of the total angular momentum J max = 
Jmax+^max is determined by the BF to SF transformation (31). The overlap and Hamiltonian 
matrices are stored on the hard disk for subsequent transformation to the SF frame. At large 
hyperradius, the couplings between different arrangements become negligibly small, and the 
Hamiltonian and sector-to-sector overlap matrices become sparse. We therefore use sparse 
matrix storage and retain all matrix elements with the absolute magnitude larger than 1CT 4 
cm -1 at p > 6.6 a Q . At each propagation step, the field-independent Hamiltonian and 
overlap matrices are constructed from the BF matrices using Eqs. (36) and (38). The 
matrix elements of the interaction with electric fields are computed directly in the SF basis 
using Eqs. (43) and added to the field-independent Hamiltonian. 

In order to avoid the overcompleteness problem 77 at small hyperradius, the eigenvalues 
of the overlap matrix smaller than the tolerance parameter e — 0.1 are discarded. The 
remaining eigenvectors of the overlap matrix are used to transform the Hamiltonian matrix 
from the primitive to orthogonalized basis (36) as described above [see Sec. II B, Eq. 
(23)]. The eigenvalue problem (25) is solved for the transformed Hamiltonian and the T ni {p) 
coefficients are used to assemble the sector-to-sector transformation matrix via Eq. (27), 
which is used to transform the log-derivative matrix to the next propagation step (54). At 
the end of the propagation, the log-derivative matrix is transformed back to the SF primitive 
basis (20) and then to the asymptotic basis (i.e. the representation in which the asymptotic 
Hamiltonian is diagonal 49 ). After these transformations, the wave function of the reaction 
complex is matched to the asymptotic Jacobi functions to yield the reactance K-matrix as 
described in the Appendix. The cross sections and probabilities for the reaction are obtained 
from the K and ^-matrices using the expressions derived in Sec. II F. 

For the LiHF system we employed the most recent potential energy surface of Aguado, 
Paniagua and Werner 78 , and used the dipole moments 6.33 D for LiF and 1.83 D for HF from 
Refs. 82 . Our largest primitive basis (20) included the vibrational states of LiF and HF with 
^LiF < 8 and whf < 3, augmented with 8 rotational states (j < 7) and 5 partial waves (£ < 4), 
which resulted in 1920 coupled channels for M — 0. The number of scattering channels is 
much larger in the SF basis, and in order to make the reactive scattering calculations feasible, 
we had to reduce the basis set parameters recommended by Week and Balakrishnan 23 for J = 
calculations. We performed the calculations in a cycle over the total angular momentum 
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projection in the range M = — 4. 

In order to test our program, we computed the cross sections for the Li + HF reaction at 
zero electric field with two different methods: (i) the SF uncoupled formalism described in 
this paper, and (ii) the standard method based on the total angular momentum representa- 
tion in the BF frame 53 ' 59 ' 65 ' 77 . The cross sections computed with method (ii) were compared 
with the calculations using the ABC code of Manolopoulos and co-workers 65 , and a good 
agreement was found for the total reaction probabilities at all J. The cross sections obtained 
with method (i) were summed over £ and the cross sections obtained with method (ii) were 
summed over J. In both calculations, we used the maximum number of rotational states 
jmax = 2 and included 6 partial waves (£ = — 5) in basis (i) and 5 total angular momenta 
(J = — 5) in basis (ii). The propagation parameters were (in units of do): p m in = 3.4, 
Pmax = 30.5, and Ap = 0.02. Figure 1 shows that the total reaction cross sections calculated 
with methods (i) and (ii) are in good agreement, which demonstrates that both the SF and 
BF formalisms are implemented correctly. Note that the SF basis is restricted to £ max = 5, 
whereas the BF cross section for J = 4 contains contributions from £ = 6. This leads to a 
small discrepancy between the SF and BF results at a collision energy of ~ 0.2 K where the 
J = 4 contribution is most significant. 

B. Numerical results 

Electric fields couple different rotational states of the reactants and products and lead 
to the formation of the field-dressed pendular states 7 . The strength of the molecule-field 
coupling can be quantified as Ed/B e , where B e is the rotational constant of the molecule. 
The rotational constant of HF (20.96 cm -1 ) is large compared to that of LiF (1.35 cm -1 ), 
and the dipole moment of HF (1.83 D) is significantly smaller than that of LiF (6.33 D). 
Therefore, the molecule-field coupling is much stronger in the entrance channel of the LiF + 
H — > Li + HF reaction. Figure 2 shows that in an electric field of 200 kV/cm, the Stark shift 
of the v — l,j — level of LiF amounts to 14.5 cm -1 , whereas the ground ro- vibrational 
state of HF is shifted only by 0.3 cm -1 . 

We consider the reaction of LiF molecules in the lowest-energy Stark state, which corre- 
lates to the state \v = 1, j = 0) in the zero-field limit. The energy of this state decreases 
with increasing the field. As the rotational levels of the HF product are only slightly mod- 
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ified by the electric field (Fig. 2), the exothermicity of the LiF(v = 0,j = 0) + H reaction 
forming HF(V = 0,f = 0) molecules decreases from 240.57 cm -1 at zero field to 226.39 
cm -1 at E — 200 kV/cm. This example shows that the exothermicity of state-resolved 
chemical reactions can be controlled with electric fields. If the energy defect between the 
initial and final ro-vibrational levels is small, it may be possible to close or open a particular 
reaction channel by varying the electric field strength. For example, the reactive channel 
\v — l,r = 2,M T = 2) — > \v' = 0,r' = 3) shown in Fig. 2 becomes closed as the field in- 
creases from zero to E = 150 kV/cm, and reopens at E > 200 kV/cm (the reader is reminded 
that the field-dressed state \tM t ) corresponds to the field-free rotational state \jMj) in the 
limit of zero electric field). As in the case of photodissociation 80 and predissociation 81 , this 
suggests that reactive scattering cross sections near threshold can be efficiently manipulated 
with electric fields. 

Figure 3 (upper panel) shows the total reaction cross section as a function of the collision 
energy. The cross section at zero field is small due to the suppression of tunneling of the 
heavy F atom under the reaction barrier 22 ' 23 . Figure 3 illustrates two important observations. 
First, electric fields enhance the reaction probability by several orders of magnitude over a 
large interval of collision energies. The largest effect is observed in the s-wave scattering 
regime, where the electric field modifies the absolute magnitude of the cross section, but not 
the dependence on the collision energy. Second, a broad resonance which appears at E c ~ 2 
cm -1 for zero electric field is completely suppressed at E = 200 kV/cm. As we demonstrated 
earlier for spin-changing collisions of CaD molecules 51 , electric fields induce the off-diagonal 
£ — > £ ± 1 transitions, which alters the relative contribution of different partial waves to the 
total cross section and suppresses shape resonances 50 ' 51 . 

An additional feature of the LiF + H reaction is the presence of rovibrational relaxation 
channels LiF(t> = l,j — 0) + H — > LiF(t>' = 0, f) + H competing with the chemical reaction. 
The lower panel of Fig. 3 shows that the vibrational relaxation is more probable than the 
chemical reaction at low collision energies. This result is consistent with the calculations 
of Week and Balakrishnan 23 . We note that because of the shape resonance at ~2 K, the 
cross section for vibrational relaxation displays the same energy dependence as the reactive 
scattering cross section for collision energies above 1 K. The coupling between different 
partial waves induced by the electric field modifies the shape resonance and suppresses both 
the inelastic and reactive cross sections in the multiple partial-wave regime. Thus, electric 
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fields can be used not only to change the absolute magnitude of the cross sections, but also 
to modify their dependence on the collision energy. We note that s-wave scattering begins to 
dominate the cross sections for vibrational relaxation at higher collision energies, especially 
at zero electric field where the reaction probabilities are very small. The presence of electric 
fields enhances the reaction probabilities in the limit of s-wave scattering, which shifts the 
upturn of the total reaction cross sections to higher collision energies. The lower panel 
of Fig. 3 shows that electric fields stimulate vibrational relaxation in the s-wave regime, 
although the effect is not as significant as for the chemical reaction. 

To elucidate the effects of electric fields on vibrational relaxation, we calculated the cross 
sections for the |i> = l,r = 0) — > \v' = 0) transition in a simpler non-reactive collision 
system CaD + He. The calculations were performed as described in Ref. 50 . Figure 4 shows 
that electric fields stimulate vibrationally inelastic scattering, in agreement with the results 
of the calculations for the LiF + H reaction (Fig. 3). The efficiency of vibrational relaxation 
is determined by the coupling of the ground rotational state of the v — 1 level with the 
rotational states in the v — manifiold induced by the anisotropy of the interaction potential. 
In the absence of an electric field, the dominant transition \v — l,j — 0) — > \v' = 0,f = 1) 
is induced by the leading anisotropic term V\ of the Legendre expansion of the interaction 
potential 

x 

Electric fields couple the initial state \v = 1, j = 0) with the rotationally excited states. 
As a result, more rotational levels within the v — 1 manifold are coupled by the V\ term 
and the coupling between the ground and the first excited vibrational states increases. To 
verify this mechanism, we calculated the cross sections for vibrational relaxation without 
the \v = l,j = 1) level in the basis set. Figure 4 shows that the cross sections calculated 
with the modified basis do not change with increasing electric field. This indicates that 
the coupling between the ground and the first excited rotational states plays a key role 
in stimulating vibrationally inelastic collisions with electric fields. The reactive scattering 
cross sections display a similar behavior (see Fig. 3), which suggests that the mechanisms 
for the electric field enhancement of chemical reactions and vibrationally inelastic collisions 
are similar. 

Figure 5 illustrates that the rate constants for the LiF + H — > HF + Li reaction increase 
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by several orders of magnitude with increasing field. The increase is not monotonic: the 
cross section for E = 150 kV/cm is one order of magnitude lager than for E = 200 kV/cm. 
This indicates the presence of an electric-field induced resonance similar to those observed 
by Avdeenkov and Bohn 83 and in our previous work 84 for non-reactive scattering. This new 
class of reactive scattering resonances is very interesting as it might allow for controlling 
chemical reactions with external fields. 

Figure 6 shows the total cross section for the Li + HF — > LiF + H reaction as a function 
of the collision energy and the electric field strength. The cross sections in the s-wave regime 
(E c ~ 1 mK) are almost unaffected by electric fields. The effects of electric fields are more 
pronounced near a scattering resonance at 0.5 K, where the cross sections are suppressed by 
a factor of 10. As discussed above, this suppression occurs as a result of the electric field- 
induced mixing of different partial waves 51 . The reaction probability at low temperatures is 
determined by tunneling under the reaction barrier 79 and the relative change of the reaction 
exothermicity induced by electric fields is small (see Fig. 2). The effect of electric fields on 
the structure of HF is relatively weak as the dipole moment of HF is quite small and the 
rotational constant of the molecule is large. We conclude that the total reaction probability 
is insensitive to the electric-field-induced interactions in the exit reaction channel. This is in 
contrast with the effects of electric fields in the entrance reaction channel discussed above. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have developed a quantum mechanical theory of reactive scattering in 
the presence of an external electric field. The approach is based on the FD hyperspherical 
coordinates and the adiabatic partitioning of the reactive scattering Hamiltonian. The total 
wave function is expanded in the eigenfunctions of the adiabatic Hamiltonian, which includes 
the molecule-field interactions. The expansion coefficients as functions of the propagation 
variable p are determined from the solution of the coupled-channel equations. The adiabatic 
eigenfunctions are constructed on a grid of p sectors by solving the adiabatic eigenvalue 
problem (13). The fully uncoupled SF basis of spherical harmonics and hyperspherical ro- 
vibrational basis functions in each chemical arrangement is used to expand the adiabatic 
eigenfunctions (Sec. II C). 

The matrix elements of the field-independent part of the adiabatic Hamiltonian (14) are 
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first computed in the total angular momentum representation using the BF basis functions 
(32). The integrals are then transformed to the SF representation using Eq. (31). The 
matrix elements of the interaction with electric fields are computed directly in the SF frame 
[Sec. II D, Eq. (40)] and added to the field-independent part (36) to yield the matrix of 
the adiabatic Hamiltonian in the SF primitive basis [Eq. (25), right-hand side]. The dipole 
moment of the atom-molecule system can be represented as a sum of the dipole moments 
of the individual diatomic molecules and a three-body term, which vanishes as the atom- 
molecule separation increases. The three-body contribution is negligible outside the strong 
potential coupling region. Inside this region, the interaction of the three-body term with 
electric fields is small compared to the atom-molecule interaction. Neglecting this term 
is therefore equivalent to ignoring the coupling between different arrangements due to the 
electric field (Sec. II D), and is a good approximation 57 . 

The matrix of the adiabatic Hamiltonian is transformed to the orthogonalized basis using 
Eq. (25) and subsequently diagonalized to yield the adiabatic eigenvalues and eigenfunctions 
at each p. During the propagation from small values of p to the asymptotic region, the wave 
function (or its logarithmic derivative) is transformed from one sector to another using the 
sector-to-sector transformation matrix (27). At the end of the propagation, the matrix of 
solutions is transformed to the asymptotic basis which diagonalizes the interaction with 
electric fields (46). The asymptotic matching procedure described in the Appendix yields 
the reactance matrix K and the scattering matrix S as well as the integral cross sections 
for reactive scattering as functions of the collision energy and the electric field strength 
(Sec. II F). 

The theory presented in this work is completely general and can be applied to any ab- 
straction reaction involving polar molecules. Because the field-independent Hamiltonian 
(14) remains the same for any atom-diatom chemical reaction 85 , the equations derived in 
Sec. II provide a general framework for including the effects of external electromagnetic 
fields in reactive scattering calculations. For example, it is straightforward to modify the 
presented formalism to describe chemical reactions in the presence of an off-resonant laser 
light 15 by changing the matrix elements of the molecule-field interaction (Sec. II E). Our 
time-independent approach can be generalized to study chemical reactions in the presence of 
microwave laser radiation 86 or radio-frequency fields 87 using the dressed-state formalism 88 ' 89 . 
Some of these generalizations are currently under development in our group. We note that 
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the theory presented here may not be applicable to insertion reactions involving the for- 
mation of long-lived intermediate complexes because the FD hyperspherical expansions are 
known to converge very slowly for this type of reactions. One possible extension of the 
present work would be to develop a formalism based on symmetric hyperspherical coor- 
dinates of Smith and Whitten 69 . This would allow for efficient numerical calculations of 
probabilities for insertion chemical reactions in the presence of external fields. 

Based on the theory developed Sec. II, we performed preliminary calculations of cross 
sections and rate constants for the LiF(i> = \,j — 0) + H and Li + HF(t> = 0, j = 0) 
chemical reactions in the presence of an external electric field. Our calculations show that 
cross sections (Fig. 3) and rate constants (Fig. 5) for chemical reactions at low temperatures 
may be sensitive to dc electric fields of less than 200 kV/cm. Our results indicate that the 
probability for the LiF + H reaction in the s-wave regime is enhanced dramatically by 
moderate electric fields in the range 100 - 200 kV/cm (Fig. 3). Low-temperature rate 
constants of abstraction reactions such as LiF + H are typically very small, and the results 
shown in Fig. 5 suggest that electric fields can be used to stimulate chemical reactions 
in cold trapped molecular ensembles. The increase of the rate with electric field is not 
monotonic, which suggests that chemical reactions at low temperatures might be affected 
by Feshbach resonances sensitive to external fields. We will explore the effects of external 
fields on scattering resonances in chemical reactions at low temperatures in future work. 

Our analysis shows that the enhancement of reaction rates at low temperatures is due to 
electric field-induced couplings between rotational states of the reactants. This suggests that 
chemical reactions at low temperatures can be selectively tuned by microwave laser fields. 
Microwave radiation couples different rotational states of the reactants, thereby inducing 
coupling between different vibrational states and reaction channels. The tunability and 
high power of microwave lasers may allow for high selectivity of control. A major thrust 
of recent experimental work has been to produce cold and dense ensembles of a variety 
of stable molecular radicals such as NH 31 ' 90 , CaH 91 , and OH 27 ' 28 . The mechanisms for the 
electric field control of chemical reactions described in this work can be readily verified in the 
experimental work with slow molecular beams 27-29 or with molecules confined in a magnetic 
or electrostatic trap 28 ' 29 ' 31 . The reaction products in external field traps could be separated 
and detected independently using the E — H gradient balance method 5 . 
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APPENDIX A: ASYMPTOTIC MATCHING 

Because the asymptotic form of the Jacobi radial functions is well known (58), it is 
convenient to re-express the FD wave function at p — > oo in the Jacobi coordinates (55). 
The asymptotic FD wave function at the end of the propagation after the transformation to 
the asymptotic basis (see Sections IID and IIF) has the form 

*= E ^P^ ,2F ^vrM T m l {p)^ a vrM T {r a ,O a ]p)Y mi (R a ), (Al) 

a,v,T,M T £,M e 

where the field-dressed FD pendular states ^> aVT M T (r a , 9 a ; p) are given by Eq. (47). Note 
that the field-dependent expansion coefficients CjM jt TM T are the same as in Eq. (57) because 
the matrix elements of the asymptotic Hamiltonian are identical in the FD and Jacobi 
coordinates (see Sec. II D). 

Since the basis functions of different arrangements are orthonormal in the limit of large 
p, we can invert Eq. (Al) to obtain 

/f f*/ 2 
dr a / dR a / dB a \ sin 2 26 a (p 5 ^)r aVT M T (r a , 6 a , p)Y e * Me (R a ), (A2) 

where \ sin 2 29 a is the angular part of the Jacobian in the FD coordinates 52 ' 53 . Substituting 
\[> from Eq. (55) yields 

df a / dR a / d6 a \ sin 2 29 a ^* omTMT {f a ,B a \ p)Y e * Mi (R a ) 

x E E p F avTM T eM e (Ra')i J a'v'T'M^(^a')Y£'M' e (Ra')- (A3) 

a'v'r'M' T e>,M' e a ' Ta ' 

As mentioned above, in the limit of large p the basis functions have zero overlap unless 
a = a'. Using the definition (57) and the orthogonality property of the coefficients Cjm^tMt, 
the integration over f a can be carried out analytically to yield 



/ 



dr a ip avTMT (r a , a ; p)VWr'M^(ra) = S TT/ S MT M^Xavr(O a ; p)^ a v'r(r a ). (A4) 
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The integration over R a is straightforward because the expansions given by Eqs. (55) and 
(Al) contain the same spherical harmonics. The matrix element in Eq. (A3) reduces after 
some algebra to the sum of one-dimensional integrals 52 ' 53 

FavrMrlMtip) = ^ ^ S aa >5 TT /5M T M^U'S Me M' i 
oc'v't'M^ e',M' e 

(rv), (A5) 

Jo 

where the renormalized FD rovibrational function a „ T (^ a ; p) is given by Eq. (30). 

Instead of Eq. (58), it is convenient to use the real boundary conditions for the Jacobi 
radial functions 

F(-Rq) = J(Ra) - N(i? a )K, (A6) 

and their derivatives 

F'(R a ) = J(R a )-N'(R a )K, (A7) 

where the primes indicate differentiation with respect to R a in arrangement a. In Eq. (A6), 
K is the reactance ^-matrix, and the diagonal matrices of incoming and outgoing waves 
J(i?) and N(i?) are composed of the modified spherical Bessel functions 52 ' 53 ' 73 . Following 
Parker and Pack 52 , we substitute Eqs. (A6) and (A7) into Eq. (A5) to obtain 

F(p) = A(p) - B(p)K, (A8) 

where 

[h(p)]avTM T eM e ,a'v'T' M^' M£ = ^aa'^Tr'^ M T M> T &U'& MiM'^^ 

r /2 

X / d9 a (f) aVT (8 a ; p)[-i(R a ')]a'v'T'M^e'M^a'v'T'(r a i), (A9) 
Jo 

The matrix B(p) has the same form with the Bessel functions J(i? Q ) substituted by N(i? a ) 
Taking the first derivative of Eq. (A5) with respect to p, we obtain 

F'(P) = ^F(p)A(p) + G(p) - H(p)K, (A10) 

where the matrix G is given by 

[^{P^avTMrtM^a'v'T'M^l'M't = ^aa'^rr'^ M T M^U'^ MiM'^ 1 ^ \ d6 a (j) aVT (6 a ; p) 

JO 

\ (AH) 



x 



d[J(i4')Wr'A^'M' d£ aV r'(?V) . a 
J7^ Ka'v'T 1 {fa 1 ) COS U a i + [J{K a >) \ a ' v ' T ' M't' M' 1 Sin V a > 

aR a i 1 ar a i 



28 



and the matrix H(p) has the identical form, with the Bessel functions J(-R Q ) substituted by 
N(i? Q ). From Eqs. (A9) and (All), we can obtain the K-matrix in the Jacobi coordinates 
directly from the log-derivative matrix 52 ' 53 ' 73 ' 74 in the asymptotic FD basis 

Y-^I1B-H IYy - ^\) A-G , (A12) 



K 



2p 



where I is the identity matrix. The ^-matrix can be obtained from the i^-matrix us- 
ing the standard Cay ley transformation 52 ' 53 . As in the conventional reactive scattering 
formalism 52 ' 53 ' 65 , the projection matrices (A9) and (All) are diagonal in all quantum num- 
bers except v. 
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Collision energy (K) 



Figure 1: Total cross sections for the Li + HF — > LiF + H reaction calculated using the SF uncoupled 
representation (circles) and the conventional BF method (full line). The SF cross sections are 
computed using Eq. (62) and summed over t = — 5. The BF cross sections are computed as 
a = ^2j(2J + l)a J and summed over J = — 4. The partial BF cross sections [(2 J + l)cr J ] are 
also presented. 
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Figure 2: Energy levels of the reactants and products of the LiF(v = 1, j = 0) + H — > HF + Li 
reaction as functions of the electric field. The initial level \v = l,j = 0} of LiF is shown by the 
dashed line. The HF states correlating to the \v = 0, j = 3) level in the zero-field limit are shown 
by the dash-dotted line in the upper panel. The zero of energy corresponds to the bottom of the 
LiHF potential well 78 at zero electric field. 
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Figure 3: Upper panel: Total cross sections for the LiF(v = l,j = 0) + H chemical reaction as 
functions of the collision energy at zero electric field (circles), E = 100 kV/cm (triangles), E = 150 
kV/cm (squares), and E = 200 kV/cm (diamonds). Lower panel: Cross sections for vibrational 
relaxation LiF(t> = 1, j = 0) + H — > LiF(v' = 0) + H summed over all final field-dressed rotational 
states as functions of the collision energy at different electric field strengths. 
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Figure 4: Cross sections for vibrational relaxation \v = l,j = Q)—> \v' = 0) in CaD + He collisions 
as functions of the applied electric field: fully converged calculations (squares), and the results 
obtained without the j = 1 level (triangles). The cross sections are calculated for M = and 
summed over all final field-dressed rotational states. The collision energy is 10~ 3 K. 
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Figure 5: Rate constants for the LiF(v = 0,j = 1) + H — > HF + Li reaction at two different 
temperatures as functions of the applied electric field. Upper panel: T = 0.05 K; Lower panel: 
T = 0.75 K. 
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Figure 6: Total integral cross sections for the Li + HF reaction as functions of the collision energy 
at zero electric field (circles) and E = 200 kV/cm (diamonds). 
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